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Abstract. Observations have revealed a large variety of structures (global asymmetries, warps, belts, rings) and 
dynamical phenomena ("falling-evaporating bodies" or FEBs, the "/3 Pic dust stream") in the disc of f3 Pictoris, 
most of which may indicate the presence of one or more planets orbiting the star. Because planets of j3 Pic have 
not been detected by observations yet, we use dynamical simulations to find "numerical evidence" for a planetary 
system. We show that already one planet at 12 AU with a mass of 2 to 5 Mj and an eccentricity 5s 0.1 can 
probably account for three major features (main warp, two inner belts, FEBs) observed in the /3 Pic disc. The 
existence of at least two additional planets at about 25 AU and 45 AU from the star seems likely. We find rather 
strong upper limits of 0.6 Mj and 0.2 Mj on the masses of those planets. The same planets could, in principle, 
also account for the outer rings observed at 500 - 800 AU. 
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1. Introduction 



Celestial mechanics - Minor planets, asteroids - Methods: N-body simulations - Stars: individual: f3 



Sinc e the discovery of the circumstellar disc of (3 Pictoris 
by ISmith fc Terrilg (1984), it beca me the most observed 



Lagrange et al. 2000, 



and best studied debris disc (see 
and references therein). However, a long-standing question 
whether (3 Pic also hosts planets, remains unanswered. 
Being an A5V star, (3 Pic is a difficult target for the ra- 
dial velocity measurements: the currently achieved preci- 
sion of hundreds m s -1 barely excludes the presence of a 
10Mj planet at 1 AU (jGalland et al.ll2006l ). The promi- 
nent edge-on disc rules out direct imaging. Transits are 
not promising either because of their low probability. At 
the same time, there is a growing bulk of in direct evidence 
for th e presence of planets in the system. IMouillet et al 



(|1997|) showed that a planet with an orbital inclination 
of 3° to 5° and a mass ranging from 0.6 Mj and 18 Mj 
between 20 and 3 AU could be responsible for the ob- 
served warp in the disc. New H ST/STIS observations of 
the warped disc bv lHeap et al. (|2000| ) chan ged these es- 
timates only slightly. lAugereau et al.l (|200lh pointed out 
that the same pla net could explain the but terfly asym- 
metry of the disc. iBeust fc Morbidellil (|2000f) found that 
a jovian planet at w 10 AU with a moderate eccentric- 
ity of e w 0.05 can explain both the warp and the ob- 
served phenomenon of falling-evaporating bodies (FEBs) . 
Krivov et al.l (|2004T ) argued that such a planet is needed 
to explain the so-called "/3 Pic dust stream" detected by 
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Table 1. Observed belt structures in the (3 Pic system 



label 


position 


ref. 


A 


w 6.4 AU 


2 


B 


« 16 AU 


1,2,4 


C 


« 32 AU 


1,2,4 


D 


« 52 AU 


1,3,4 



Ref erences: 

l= flWahhai et al.ll2003h , 2= flOkamoto et al.ll2003) . 

3= llTelesco et al.ll2005l ). 4= l!Golimowski et aril2006f) 



meteor radar AMOR ( Baggalev 2000f) . Finally, most re- 
cent observations revealed several belt-like structures in 
the inner disc (Table [l| that can be attributed to the 
presence of planets. This work is an attempt to constrain 
the parameters of suspected planets in the light of these 
observations. 

2. Presumed planetesimal belts 

Okamoto et al. I (|2004h performed high-iesolution spectro- 
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scopic observations in the 10-/mi band to identify con- 
centrations of submicron-sized silicate dust at 6.4 AU, 16 
AU and 30 AU and interpreted these by dust-replenishing 
planetesimal belts at those locations. To check this inter- 
pretation, we have made test runs of our collisional code 
(Kri vov et al. 2006). The code enables simulations of a 
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circumstellar disc of solids over a wide range of sizes - 
from planctcsimals to fine dust - taking into account stel- 
lar gravity, radiation pressure, as well as destructive and 
cratering collisions. We took a planetesimal belt of ob- 
jects with radii from 0.15/mi to 7 km and a total mass of 
O.3M0, semimajor axes from 6 and 7 AU and eccentricities 
between 0.0 and 0.1 and evolved it to a quasi-steady state 
to obtain a spatial distribution of dust material sustained 
by the belt. We then calculated the blackbody thermal 
emission of that dust and took a standard line-of-sight in- 
tegral to obtain a brightness profile as a function of the 
projected distance from the star. The mass an d luminosity 



of Pic were taken to be 1.75M Q and 8L Q IjCrifo et al 



19971) . We assumed a bulk density of solids of 2.5 g cm . 

The results are shown in Fig. [TJ The flux is dominated 
by emission of grains in bound orbits, or a-meteoroids. 
However, particles in hyperbolic orbits, or /3-meteoroids, 
which are particles smaller than w 2/im, also make a sen- 
sible contribution. The "half-peak" brightness is about 
a factor of two larger than in the adjacent parts of the 
ring, roughly consist ent with the observations (Fig. 2 of 
Okamoto et al.ll2004j ). A more accurate comparison is not 
possible, for both a large difference between the 6.4 AU 
SW and NE peaks and the fact that the peak brightness 
depends on the assumed width of the planetesimal ring. 
The total mass of the planetesimal belt of 0.3M ffi , compa- 
rable to that of the Kuiper belt in the solar system, leads 
to the peak brightne ss ~ 0.3 Jy arcse c" 2 , which is close to 



the observed values (jOkamoto et al.ll2004T ). Therefore, we 



can conclude that the observed brightness peak at 6.4 AU 
can indeed be attributed to an invisible planetesimal belt 
at approximately the same location. 

The results for the outer belts look similar and there- 
fore are not shown here. Furthermore, our collisional mod- 
elling implies that, with a reasonable accuracy, possible 
interaction between the belts can be neglected, and that 
they can be treated separately. 



3. Presumed planets 

3.1. Numerical model of the planetesimal disc 

Numerical simulations of the motion of planetesimals were 
carried ou t with the mercury6 integration package by 
Chambers! (Il999h . The planetesimal disc was modeled by 
a set of massless particles, initially distributed in equidis- 
tant, circular and plain orbits around the star, in which 
we placed one or more planets. Since we want to focus our 
study on the inner part of the disc, the orbits of particles 
ranged from 1 AU to 70 AU with a step size of 0.1 AU. All 
particles that were in trojan-type motion with one of the 
planets were removed. Altogether, we performed about 
100 different simulations and in each run the disc con- 
tained 20730 particles. The integration interval was set to 
12 Myr in accordance w i th the supposed ag e of the system 
(jZuckerman et alj|2001i lOrtega et al.ll2002n . 
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Fig. 1. The 10/xm edge-on brightness of a dust "subdisc" 
produced by a planetesimal belt between 6 and 7 AU (as 
shown with vertical lines). Dashed and dotted lines: con- 
tributions from a— and /3-meteoroids, respectively; solid: 
their sum. 



3.2. One planet 

The first s imulations were carrie d out with one planetary 
perturber. lOkamoto et al.l (|2004h suggested that the bor- 
ders of the belt at 6.4 AU are created by resonances with a 
planet at 12 AU, similar to the main belt in our solar sys- 
tem which is confined by a 1:2 and 1:4 mean-motion reso- 
nance (MMR) with J upiter. Thus we adopted a = 12 AU. 
Mouillet et alJ (jl997l ) give a possible range for the incli- 
nation of a planet between 3° and 5° so that we chose 
i = 4°. Three different values of the eccentricity were 
tested: e = 0.01, e = 0.1 and e = 0.2. We tried the planet 
with several masses between 2 and 5 Mj, the values which 
roughly bracke t the range sugge sted to explain the warp 
(see Table 4 of lHeap et ai1l200ul) . 

The simulation results are shown in Fig. [2j Note that 
this and other plots are meant to show the positions of 
belts and gaps only; to determine their strength one would 
need to convolve the distributions with the assumed ra- 
dial density profile of the disc. The planet with m = 2Mj 
opens a gap in the disc and causes a peak at 6.4 AU and a 
smaller peak close to 16 AU, corresponding to the A- and 
B-belts. The eccentricity has no influence on the peak at 
6.4 AU but the peak at 16 AU shifts outwards to « 20 
AU for e = 0.2. For the planet with m = 5Mj, the peak 
at 6.4 AU is also present, but the second peak is at about 
20 AU rather than 16 AU for all three values of eccentric- 
ity. The error bars for the estimated parameters of this 
planet are listed in Table H The limits were determined 
by changing the semimajor axis, eccentricity and mass and 
by checking, whether the resulting peaks in the distribu- 
tion of test particles still satisfy the reported observations. 
To avoid possible confusion in interpreting the uncertain- 
ties of semimajor axes of the proposed planets, we note 
that we assumed here the positions of the peaks found 
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Fig. 2. Effect of a planet at 12 AU with m = 2Mj on the 
inner part of the disc of (3 Pic (solid line: e = 0.01; dashed 
line: e = 0.1; dashed-dottcd line: e = 0.2). 
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Fig. 3. Time needed to reach 1000 AU by the scattered 
particles. 



by lOkamoto et al.l (|2004) to be exact and give in Table [2] 
and other places of the paper the estimated intri n sic er - 
ror of our simulations. However, lOkamoto et al.l (|2004f ) 



have sampled the disc with a step of 3.2 AU, performing 
their measurements at discrete distances from the star of 
3.2 AU, 6.4 AU, 9.5 AU, and so on. Therefore, the actual 
uncertainty of the semimajor axes of the predicted plan- 
ets is determined by the uncertainty of the peak locations, 
re 2 to 3 AU. 

Although we propose that a second planet around (3 
Pic makes the peak near 16 AU more prominent (see 
Sect. I3.3p . already this planet alone could produce a no- 
table peak. There is an analogy with the Hilda group 
of asteroids in our solar system that move near the 2:3 
resonance with Ju piter and Plutin os in 3:2 MMR with 
Neptune (see, e.g.. lMorbideii3 12002). In our case, the 3:2 
MMR with the planet at 12 AU is located at ~ 15.7 AU. 
Many planetesimals in this region are resonant-protected 
from close encounters with the planet and form a peak 
with respect to the non-resonant background population. 
Similarly, the inner peak at 6.4 AU is strongly "supported" 
by the 2:5 MMR in much the same way as the Koronis as- 
teroid family by its 2:5 MMR with Jupiter. 

The particles originally placed near the planet were 
scattered out of their position; some of them were thrown 
out of the system. Fig. [3] shows the time such a particle 
needs to reach a distance of 10 3 AU. The majority need 
~ 0.5 Myr to reach this distance and thus they might form 
some observable features in the outer parts of disc, for 
insta nce the rings observed at 500 - 800 AU (jKalas et al 



20011 ). Discrete rings may be a result of ejection of several 
large planetesimals. Fig.|4]shows, for three of the particles, 
how close encounters with the planet at 12 AU can change 
successively the apastron distance until it lies in the outer 
part of the disc. 



5 1 200 p 

V 1000 - 

| 800 - 
| 600 - 
§ 400 - 



« 180OOO 



200000 



300000 400000 
time [y] 



500000 



600000 



Fig. 4. Evolution of the apastron distance of three scat- 
tered particles 



Beust k Morbidel li (2000) analyzed mechanisms pro- 
posed to explain the observed phenomenon of falling- 
evap orating bodies (FE Bs). All of these - c lose encoun- 
ters ( Beust et al. I ll99ll) . Ko zai mechanism dBailev et al 



I992I) , trapping in MMRs (|Beust k Morbidellil ll996) 



involve the presence of at least one large planet. They 
showed that a planet with ps 2 Mj at pa 10 AU with a low 
eccentricity can account for the detected infall of small 
bodies onto the star. Furthermore, they noticed that the 
parameters of the pla net lie well inside the limits given by 
iMouillet et al. (1997) for a planet that can cause the ob- 
served warp. Because these parameters are close to those 
we found above, we can conclude that already a single 
planet with a mass of m w 2Mj at 12 AU, e ^ 0.1 and 
i ps 4° can probably account for three of the major dynam- 
ical phenomena observed in the disc of (3 Pic: the warp, 
the A- and B-belts and the FEBs. 



3.3. Two and three planets 

The second planet can be invoked to account for the peak 
around 32 AU (C-belt) that could not be created by the 
planet at 12 AU. If we follow the same strategy as for 
the first planet, we find that the best fit is a planet with 
m = 0.5Mj at a = 25 AU and e = 0.01 that clears a 
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Fig. 5. Influence of two planets at 12 and 25 AU with 
m = 2M,j and m = 0.5Mj (both with e = 0.01) on the 
disc. 
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Fig. 6. Effect of three planets at 12, 25 and 44 AU with 
m = 2Mj, 0.5Mj and O.lMj (all with e = 0.01) on the 
disc. 



gap and causes a new peak at 32 AU. The distribution of 
particles in the two-planet case is shown in Fig. [5] 

Similarly, a third planet is needed to explain the peak 
around 52 AU (D-belt). The best fit is with m = 0.1 Mj at 
44 AU. Figure [5] shows the result of our simulations with 
three planets. The D-belt is now created, and the peak at 
32 AU became more prominent because it is now bordered 
by a planet on both sides. Also, the peaks at 6.4 and 16 AU 
are still present. It is important to mention that a more 
massive second planet, with a mass larger than rs 0.6Mj, 
can be excluded because it would destroy the belts. The 
same applies to the third planet which cannot be more 
massive than w 0.2Mj. Table [2] summarizes the best-fit 
parameters of all three planets and their errors bars. 



Table 2. 

system 



Parameters for the proposed (3 Pic planetary 



Planet 


m [Mj] 


a [AU] 


e 


1 


z - u -0.5 


12 ± 0.5 


01+ 01 

u - ui -0.01 


2 


0.5 ±0.1 


25 ± 1 


o.oil8:8S 


3 


1+ 01 

U - X -0.03 


44 ± 1 


o.oil8:8? 



4. (3 Pic - a resonant system? 

The orbital periods of the proposed three planets are close 
to rational commensurabilities. Namely, the periods of 
planet 2 and planet 1 are nearly in a 3:1 ratio; those of 
planet 3 and planet 1 are close to a 7:1 ratio; and planet 3 
and planet 2 are near a 7:3 commensurability. Planets in 
some of the extrasolar planetary systems discovered so far 
are known to be locked in mean-motion resonances (e.g. 
GJ876 or 55 Cnc). To see if the possible planets of (3 Pic 
really show resonant motion, we have chosen the strongest 
of the three resonances mentioned above, the 3:1 MMR be- 
tween the two inner planets, and calculated the resonant 
angle (|ji et al.ll2003h 



= Ai - 3A 2 + (cJi + cJ 2 ) , 



(1) 



where A^ and c5,; are the true longitude and the longitude 
of periastron of the i-th planet. 

Figure[7] (top) shows the case, where planet 2 is located 
at the exact position of the resonance (a± — 12.00AU and 
a 2 = 24.96AU); both planets have an initial eccentricity 
of 0.1. Other parameters are the same as in the simula- 
tions described in Sec. 13.31 The resonant angle librates 
around 0°, indicating the resonant locking. However, the 
eccentricties are not strongly affected by the resonance 
which seems to be a protective one, as it is the case for 
GJ876 and 55 Cnc. That the effect is weak is not surpriz- 
ing: the resonance is rather shallow (as seen from a rather 
large libration amplitude), which, in turn, traces back to 
moderate masses of both planets and a large separation 
between their orbits. 

We have found that a decrease in semimajor axis of the 
second planet by sa 0.1 AU or an increase by « 0.2 AU is 
enough to transfer 9 from libration to a circulation mode 
(Fig. bottom) . This value can be compared to analytic 
estimates of the resonance width. The semimajor axis of a 
body locked in a resonance oscillates between a — Ao max 
and o + Aomay, where the libration amplitude is given by 



(jMurrav fc Dermottl[l999l ) 



"2 



16 mi 

3 m* 



P 



p + q 



fae-2 



with / d = 0.5988 for the (p+q) : p = 3 
Assuming e 2 = 0.1, this leads to 



(2) 



1 outer resonance. 



0.3AU, 



(3) 
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Fig. 7. Resonant angle 9 for the 3:1 resonance of planet 1 and planet 2 (left) and evolution of the eccentricities (right). 
Top: the case were planet 2 is located at the exact position of the resonance. Bottom: planet 2 is shifted by 0.1 AU 
inwards. Both planets have an initial eccentricity of 0.1 and an inclination of 4°. 



which is close to the value we found empirically. Thus the 
resonance width is smaller than the uncertainties in the 
semimajor axes of the planets (Aa ~ 1 AU for our fits ; 
Aa ~ 3AU for the original observations of Okamoto et al 



(2004)). 



The eccentricity plays an important role as well. Our 
best fits for the three planets imply small initial eccen- 
tricities (e ~ 0.01). With these small values, the resonant 
argument in our simulations was always circulating rather 
than librating. Libration of the resonant angle was ob- 
served, starting from e *J 0.07. The fact that the reso- 
nance gets thinner at low eccentricities is also seen from 
Eq. [21 Aa max oc -y/ej. Since the uncertainties of the de- 
rived semimajor axes are much larger than the precision 
required to distinguish resonant from non-resonant config- 
urations, we cannot conclude whether the possible (3 Pic 
system is a resonant one or not. As mentioned above, some 
exoplanetary systems are known to be resonant. It is pos- 
sible that a planetary system emerges as a resonant sys- 
tem — for instance, as a result of differential migration in 
a gas eous disc (e.g., iFerraz-Mello et aL 20051 : Kiev et al. 
l2005h . Alternatively, the /3 Pic system could resemble the 
Jupiter-Saturn configuration in our Solar System. Here 



the two planets are close to, but not locked in, a 5 : 2 
resonance. 



5. Conclusions 

We have used numerical simulations to investigate the ef- 
fect of one or more planets on the disc of (3 Pic. The goal 
was to find a minimum set of planetary perturbers which 
could be responsible for as many features/phenomena ob- 
served in the disc as possible. 

One important result is that already one planet at w 
12 AU with a mass of m ss 2Mj and an eccentricity of 
e ^ 0.1 is able to cause three major known features: (i) 
two of the four belt-like structures listed in Table [TJ (ii) 
th e parameters of t he pla net lie well inside the limits given 
by iMouillet et al. ( 1997 ) for a planet responsible to the 



wa rp; and (iii) this pla n et is similar to the one proposed 



by iBeust fc Morbidellil (|2000h to explain the FEB phe- 
nomenon (a ps 10 AU, m = 2Mj, low e), so that the two 
can be considered identical. Another result is that two 
additional, more distant planets would naturally explain 
further planetesimal belts suggested by some observations. 
We find rather strong upper limits on the masses of those 
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planets: s» 2Msaturn for a second planet at 25 AU and 
ph 4Mneptune for a third one at 44 AU. More massive 
perturbers would destroy the belts. Figure [5] shows the 
positions of the proposed planets in the semimajor axis- 
mass plane. Overplott ed are regions exclud ed by radial ve- 
locity measurements ( Galland et al. 20061 ). as well as the 
phase-space locations of the warp-inducing planet and the 
"FEB planet" . Interestingly, all three planets are close to 
the "warp planet line" . 

The same planet(s) could, in principle, also account for 
the outer rings ob served at 500 - 800 AU from the star 
( Kalas et al.ll200lt ). A few large planetesimals could have 
encountered one of the inner planets in the past and, per- 
haps after tidal disruption during these encounters, been 
sent by the planets into escaping orbits or those with apoc- 
enters at hundreds of AU where the rings are observed. 
Further work is also needed to answer the question if the 
planetary system of (3 Pic is a resonant one or not. If it re- 
ally exists and the planets are in a resonant configuration, 
this may indicate that migration played an important role 
in its formation. 
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